Scalable scheme for entangling multiple ququarts using linear optical elements 
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I. INTRODUCTION 
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(Dated: to appear in Phy. Lett. A (2007), doi:10.1016/j.physleta.2007.05.033) 

We report a scalable linear optical scheme for generating entangled states of multiple ququarts 
in which the individual single-ququart state is prepared with the biphoton polarization state of 
frequency-nondegenerate spontaneous parametric down-conversion. The output state is calculated 
with the full consideration of the higher order effect (double-pair events) of spontaneous parametric 
down-conversion. Scalability to multiple-ququart entanglement is demonstrated with examples: 
linear optical entanglement of three and four individual biphoton ququarts. 

t~-l« ". PACS numbers: 03.67.Mn, 42.65.Lm, 03.67.-a, 03.67.Hk 
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In quantum information, the quantum bit (qubit) is often considered as the basic unit of information processing. 
C*~) Physically, the qubit is equivalent a two-dimensional quantum state, such as, electron spin, photon polarization, 
two discrete energy levels, etc. Due to the intrinsic quantum nature of the qubit, it has two distinct properties in 
comparison to the classical bit: quantum superposition (for single qubit) and quantum entanglement (for multiple 
O i' qubits). In quantum information, these distinct quantum properties of the qubit are harnessed to perform certain 
, 1 * , information processing and communication tasks better than using classical means. 

The dimension of the quantum system, however, need not be fixed at two. Continuous variable (i.e., infinite di- 
mensional) quantum states have long been studied in quantum optics to understand squeezing and certain quantum 
information protocols, such as, quantum teleportation and quantum cryptography have been implemented with con- 
tinuous variables [3, [2, H[ • It is then natural to consider a high-dimensional discrete quantum state as the potential 
resource or the basic information unit in quantum information. In fact, recently, there has been much interest in the 
.D-dimensional (D > 2) quantum state (quDit) for potential quantum information applications and for fundamental 
tests of the quantum theory 0, H H, 0, II • 

In the quantum optical approach to quantum information, the qubit is often implemented with a certain easy- 
to-control two-dimensional degree of freedom of a photon, such as, polarization, path, time-of-arrival, and phase 
O 1 SOS El- Similarly, the quDit may be implemented with a high-dimensional degree of freedom of a photon, such as, 
multiple paths, multiple time-bin, angular momentum states, etc. Of particular imp ortance and interest, among the 
photonic quDit implementation schemes, are reported in Ref. [l^, EH and Ref. jl4 fill]. m which the qutrit and the 
^ ququart states, respectively, are implemented with a pair of photons, rather than with a single photon. 

In the recent years, there have been a number of experimental demonstrations of two-quDit entangled states 
using the hi gh- dimensionality of photonic angular momentum, transverse momentum-position, time-of-arrival, etc. 
[HI, TitI - [H, Tl9L [20j . What is common in these seemingly different schemes is that they all utilize the inherent 
position-momentum or angular-momentum correlations exist between the two photons of spontaneous parametric 
down-conversion. As a consequence, two entangled quDits are always generated as a pair and the entanglement 
schemes are, in principle, not scalable to multiple quDits. 

In experimental quantum information research, it is desirable to be able to prepare individual quDit states separately 
and to have a scheme to entangle such individual quDits in a scalable manner, when needed. This bottom-up approach 
to entanglement generation is more or less impossible for the quDits reported in Ref. 0, [13, [H, [3, H(| due to the 
physical nature of the photonic degrees of freedom used to implement the quDit. 

Recently, we proposed a linear optical scheme for generating entangled states of two photonic ququarts in which 
the individual ququart state is made of a pair of photons of different frequencies (2l| . In this paper, we investigate 
realistic and experimentally realizable two-ququart and multi-ququart entangling schemes in which the single ququart 
state is based on frequency-nondegenerate two-photon polarization state of spontaneous parametric down-conversion 
(SPDC) [3, UH. The SPDC process generates a coherent superposition of the vacuum, the two-photon amplitude, the 
four-photon amplitude, etc [22j. Here, we discuss experimentally relevant (postselected or non-postselected) output 
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FIG. 1: The single biphoton ququart detection scheme [l4], EH- The dichroic beam splitter (DBS) splits the two incoming 
frequency-nondegenerate photons (Ai and A2) into two spatial modes. Polarization projection measurement is performed using 
a quarter-wave plate (QWP), a half- wave plate (HWP), and a polarizer (POL) on each of the photon. CC is the coincidence 
counter. 

states of the linear optical entangling scheme by considering the real SPDC photon states (i.e., coherent superposition 
of the vacuum, the two-photon, the four-photon, etc. terms) as the input [23|. We also discuss, with examples, how 
the beam splitter-based two-ququart entanglement scheme can be scaled up to generate entangled states of multiple 
biphoton ququarts. 

II. TWO-PHOTON QUQUART 

Unlike other schemes that implement quDit states on the multi-dimensional degree of freedom of a single-photon 
state 0, E3 EH [H, H3| , the ququart, a four- level quantum state, discussed in this paper is based on the two-photon 
polarization state of frequency-nondegenerate SPDC [24j. 

Enco ding multi-dimensional quantum states on a photon pair, rather than on a single photon, was first reported 
in Ref. [U 0) EH EH and there are a number of advantages to this approach. First, since the two-fold coincidence 
detection is the basic measurement for the pair photon, the scheme is less susceptible to environmental optical noises 
and the detector dark counts. Second, the photon pair state can be prepared reliably and efficiently with the SPDC 
process, while there is no efficient and reliable source of single-photon states at the moment (25|. Third, this two- 
photon approach is particularly efficient and useful for the ququart state implementation. The degree of polarization of 
frequency-nondegenerate two-photon polarization state of SPDC is not invariant under SU(2) transformations [26l |27]| . 
Therefore, all four biphoton ququart basis states can be easily accessed by a set of waveplates, i.e., linear optically 
[iH EH ■ Finally, as we discussed earlier in Ref. [2l| , it is possible to generate entangled states of two individual 
ququarts using beam splitters, such as, an ordinary 50/50 beam splitter (BS), a polarization beam splitter (PBS), 
and a dichroic beam splitter (DBS). In addition, as we shall show in section ITVl it is possible to generate entangled 
states of multiple ququarts linear optically. 

Let us start by defining the basis states that form the biphoton ququart [li], EH • For the frequency-nondegenerate 
two-photon polarization state of SPDC, it is not difficult to realize that the two-photon polarization state must belong 
to the four-dimensional Hilbert space and that the four orthonormal basis states can be written as, 

|-ffAi,#A 2 ) 

\H\ tl V\ 2 ) 
\V Xl ,H X2 ) 

where subscripts Ai and A2 label two different wavelengths of the photons. H and V refer to the horizontal and 
vertical polarization states, respectively. It is important to note that these two photons must occupy the same spatial 
mode. 

The most general form of the biphoton single-ququart state is then given by the coherent superposition of the 
two-photon amplitudes (the ququart basis states) defined in eq. (JTJ), 

|V>}=co|0)+ Cl |l)+ C2 |2)+c 3 |3), (2) 

where Cj — \cj\e lVj are the complex probability amplitudes and satisfy the normalization condition YIj—q \ c j\ 2 = !• 

The state in eq. © can be experimentally implemented using SPDC. For example, two type-I SPDC sources 
can each be used to prepare the ququart basis states \H\ 1 ,H\ 2 ) and \V\ 1 , V\ 2 ) and two type-II SPDC sources for 
each of the basis states |-ff\ij V\ 2 ) and \V\ 1 ,Hx 2 ). Coherent combinations of these four SPDC amplitudes will then 
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FIG. 2: The proposed linear optical two-ququart entanglement scheme. Waveplates WP1 and WP2 transform the biphoton 
states to the desired ququart states \4>)a and \ip)b- Three different types of beam splitters are investigated: an ordinary 50/50 
beam splitter (BS), a polarizing beam splitter (PBS), and a dichroic beam splitter (DBS). 

make up the general single-ququart state in eq. ([2]) and in this way we have the full control of the coefficients in 
eq. (J2j) . Experimentally, however, this is a bit inconvenient and complicated as coherently combining (i.e., without 
any distinguishing temporal and spectral information) type-I and type-II SPDC is certainly not an easy problem due to 
the group velocity mismatches and the two-photon bandwidth differences between the two phase matching conditions 
of SPDC. A simpler procedure would be to generate one basis state in eq. ([1} using either type-1 or type-11 SPDC 
and then to transform the state to a more complex single-ququart superposition state by using a set of waveplates, 
as reported in Ref . [TJ, . 

To detect the biphoton ququart state, the biphoton ququart detection scheme shown in Fig. Q] can be used [Til [TBI]. 
As we can see from eq. ([TJ, the two photons have the different frequencies Ai and A2 so that they can be easily 
separated into two spatial modes by using a dichroic beam splitter. The polarization measurement is performed 
using a quarter-wave plate, a half-wave plate, and a polarizer on each of the now separated single-photons and the 
coincidence events between the two detectors are recorded at the coincidence counter to obtain the joint polarization 
measurement output. To completely characterize the single biphoton ququart state (or the density matrix), total of 
sixteen joint polarization measurements should be performed [Til . [HI l2q | . 



III. ENTANGLING SCHEMES FOR TWO BIPHOTON QUQUARTS 



To discuss the linear optical two-ququart entangling scheme, let us start with two independent biphoton ququarts 
\^p) a and |?/>)b of the form, 

H) a = c \0) a + ci|l) a + c 2 |2) Q + c 3 |3) a , 

|V)6 = ^|0) 6 + ci|l) 6 + 4|2) b + 4|3) 6 . [ ' 

In the linear optical two-ququart entanglement scheme, the entangled state of two biphoton ququarts is generated 
by interfering them at a beam splitter. The schematic of the proposed linear optical two-ququart entanglement 
experiment is shown in Fig. [5] A femtosecond uv laser pulse is split into two and pump the SPDC generating BBO 
crystals. Using the waveplates WPI and WP2, the biphoton state is transformed to the desired ququart states \ip) a 
and \ip)b and they are made to overlap (or interfere) at the beam splitter. Here we consider three different types of 
beam splitters, an ordinary 50/50 beam splitter (BS), a polarizing beam splitter (PBS), and a dichroic beam splitter 
(DBS). Since Shih- Alley and Hong-Ou-Mandel type quantum interference is used to generate entangled states of the 
two ququarts, it is necessary that the SPDC process is pumped by a ultrafast laser [2J, [29|, [3(| • The delay r adjusts 
the overlap of the two ququarts at the beam splitter, hence controlling temporal distinguishability. 

The biphoton state \ip) ou t at the output modes c and d of the beam splitter is then investigated for three different 
types of beam splitters, an ordinary 50/50 beam splitter (BS), a polarizing beam splitter (PBS), and a dichroic beam 
splitter (DBS). Assuming that the beam splitter is loss-less, we can describe it as a unitary matrix acting on incoming 
photons. Therefore, once the beam splitter unitary matrix is fully determined, the output state \ip) ou t can be readily 
calculated. 

Note that we are only interested in the biphoton states that can be expressed in the ququart basis states in eq. |T]). 
This is justified because, to detect a single biphoton ququart, it is necessary to perform the two-photon coincidence 
measurement shown in Fig. [T] Any biphoton amplitudes, at the output modes of the beam splitter, that cannot be 
expressed in the biphoton ququart basis states in eq. ([T]), therefore, do not result in detection events and can be 
ignored. 
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The input two-ququart state to the beam splitter, as shown in Fig. [2j is obviously a product states. So, at first 
sight, one may simply write the input state to the beam splitter as, 

\i>)in = \i>)a <8> \lp)b- (4) 

Equation which assumes that each SPCD source emits exactly one pair of photons, however, is incorrect in 
general for the interference experiment involving two separate SPDC sources. The SPDC process, in reality, leads 
to coherent superposition of the vacuum, two-photon (single pair), four-photon (double-pair), and higher-order terms 
[13, H3|- The individual ququart states entering the BS, therefore, should be modified as, 

\^)SPDC _ | woc ^ + + ^2^ + . . . ? 

where |t>ac) is the vacuum, 77 is related to the two-photon (ququart) event efficiency, and \ip) a (\ip)b) is as defined 
in eq. ([T]). Note that the first order terms of rj are the desired two-photon event terms. In two-photon coincidence 
experiments, higher order terms {rj 2 and higher) are normally dropped since rj <C 1. 
Clearly, the input state to the beam splitter in eq. (|3|) should also be modified as 

wi PDC = w a pDc ® m pDC . (6) 

Substituting eq. ([5]) into eq. ©, we find that, 

where \ip) a \ip) a , for example, is a short-hand notation for \ip) a (g) |^)„ (g) |wac)b which refers to the event that four- 
photons are generated at the SPDC source a while no photons are generated at the SPDC source b. Here, we have 
dropped the non-detectable vacuum component \vac) as well as amplitudes that are higher order than r/ 2 . 

Note now that we can no longer ignore the four-photon event of the SPDC, the second and the third terms in 
eq. (J7J), as they have the same probability as the desired input state, the first term of eq. ([7]). 



A. Ordinary 50/50 beam splitter 

For an ordinary 50/50 beam splitter (BS), the unitary transformation is obtained from the welll- known beam splitter 
input-output relation, 

Cjk = (a jk +ib jk )/V2, 
djk = (ia 3 k + b jk )/v2, 

where c, for example, is the annihilation operator for a photon in mode c. The subscripts j and k refer to the 
polarization mode (\H) or |V)) and the frequency mode (Ai or A2) of the photon, respectively. 

For a single-photon at the input mode of the BS, as in Fig. [2 there are four possible states (two spatial modes a 
or b and two polarization modes H or V). We can therefore define the single-photon states as follows [2l| 

a] ?fe |0) = (1,0,0,0) T , a^|0) = (0,1,0, of, 
S^ fe |0) = (0,0,l,0) T , 6|, fc |0) = (0,0,0, if. 

Equation ([9|) then forms the basis set for the unitary transformation U^ 3 of the single-photon of frequency Ai due 
to the BS. Using eq. JHJ and eq. ([9]), we find that U^f is given as, 



U 



BS 
Ai 



1 



/ 1 i 
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1 1 
\0 i 1 



(10) 



Note that U^ 3 = U^ 2 S as the BS unitary transformation does not depend on frequency and polarization states. 

Remember now that the ququart is in fact made of a pair of photons. The two-ququart interference-entangling 
scheme shown in Fig.(2J therefore, can be understood as a four-photon interference scheme with four photons incident 
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at the input ports a and b of the BS, according to eq. (JT]). The unitary transformation that describes this effect is 
then calculated to be, 

U BS = U^ s <8> C/f/ ® C/f/ C/f/. (11) 

Equation (fTTj) , a 256 x 256 matrix, describes all possible unitary transformations by the BS due to the input state 
given in eq. (f4|. 

Let us first consider the state transformation of the two-photon term |i/ , )o|V')6 m ec L- ©■ Using eq. (fTT) . the output 
state of the BS for the state \ r ^) a \4')b is calculated to be, 

\lp)a\lp)b ^ ^(CQC' 3 + C 3 Cq - dc' 2 - C 2 ci) 

x (|0) c |3} rf + |3> c |0) rf - |l) c |2) d - |2) c |l) d ) (I 2 ) 

In equation (|12[) . there are total of 64 additional four- photon amplitudes that cannot be expressed in the ququart basis 
states defined in eq. ([1]). Since these amplitudes cannot be detected by the ququart detection scheme shown in Fig.[TJ 
they have been dropped out. The other 176 amplitudes simply vanishes due to the Hong-Ou-Mandel type destructive 
quantum interference. (Note that, since the state is pure, the output state is described by a 256 x 1 column vector.) 

The output states resulting from the four-photon terms, the second and the third terms in eq. J7]), can be evaluated 
in a similar way. After a lengthy calculation, we find that the output state for the second term in eq. ([7]) is given as, 

\ip)a\ip)a\vac)b - co|0) c |0) d - c Ci|0) c |l) rf - c c 2 \0) c \2) d - i(c c 3 + c 1 c 2 )j0) c |3) d 

- coCi|l) c |0)d - cf |l) c jl)d - ^(c c 3 + cic 2 )|l) c |2)d - cic 3 |l) c |3) d 

- coc 2 |2) c |0) d - i(c c 3 + Cl c 2 )|2) c |l) d - c 2 2 \2) c \2) d - c 2 c 3 |2) c |3) d ( 13 ) 

- ^( c o c 3 + cic 2 )|3) c |0)d - cic 3 |3) c |l) d - c 2 c 3 |3) c |2) d - c§|3) c |3) d 

Here, similarly to eq. (p~2|) . only the biphoton terms that can be expressed in the ququart basis state defined in eq. ([I]) 
are shown. Total of 192 biphoton amplitudes are omitted (postselected out) in eq. (fl3|) as they do not trigger the 
four-fold coincidence circuit which is at the heart of the two-ququart detection scheme. 

The third term in eq. ([7|), has the identical output state as in eq. (fT3"|) except that all the coefficients are now 
primed, following the definition of \ip)b in eq. (J3]). We shall define this output state as |V , )i 35 ■ 

The final output state of the BS postselected by the four- fold coincidence measurement, due to the input state in 
eq. j7]) , should then be expressed as the superposition of eq. (fT2|) and eq. (fT3|) , 

Clearly, the output state is a rather complicated two-ququart entangled state. Experimentally, one may want to 
generate a more manageable two-ququart entangled state and, to do so, it is necessary to make appropriate choices 
of the coefficients in eq. <j3j> - 

As an example of the two-ququart entangled state, consider the case where cq — 1 and c 3 = 1 in eq. (|3|). Then the 
real input state to the BS is given from eq. (O as 

l^>f„ PDC = |0)„|3) 6 + |0) o |0) o + |3) b |3) b . (15) 

The second and the third terms come from the double-pair events of SPDC. Given the above state as the input to 
the BS, the postselected output state is then found to be 

IV'W = -|0) c |0) d - |3) c |3) d + j(|0) c |3) d 
+ |3) c |0) d -|l) c |2) d -|2) c |l) d ), 



which is a rather complicated-looking two-ququart entangled state. 
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B. Polarizing beam splitter 

To investigate two-ququart entangling properties of the polarizing beam splitter (PBS), we first need to determine 
the proper form of the PBS transformation. Following the spatial mode structure in Fig. [2J we have, 

o-Hk — > cjjfe, avk — > dvk, , . 

(17) 

bnk —> d-Hk, bvk — > Cyfe, 

where the subscript k refer to the frequency mode (Ai and A2) of the photon. 

The PBS unitary transformation for the single-photon of frequency Ai is then given as, 



jjPBS 



/l 000 

1 

10 

\0 1 



(18) 



Since the PBS unitary transformation does not depend on the frequency mode, U BBS = U BBS . 

The PBS unitary transformation that corresponds to the input state in eq. ([7]) is therefore calculated to be, 

U PBS = U BBS ® U BBS ® U BBS ® U BBS , (19) 

which is again a 256 x 256 matrix. 

Applying the PBS unitary transformation in eq. (fT9l) to the the two-photon term in eq. ([7]), we get the following 
postselected two-ququart state (out of total of 16 amplitudes) in the output modes c and d of the PBS. 

WaWb ^ CoC |0) c |0) d + C lC ' x \l) c \l) d (20) 

+ c 2 4|2) c |2) d + c 3 4|3) c |3) d 

For the four-photon amplitudes in eq. 0, the postselected output states (out of total of 16 amplitudes) are given 

as, 

\ip)a\ip)a\vac) b PBS > 2(c c 3 + c x c 2 ) |0) c \ 3) d , 

DDC V / 

\vac) a \if>) b \tl>) b ^ 2(44 + c;4)|3) c |0) d . 

Finally, we find that, given the input two-ququart product state in eq. ijT]). the PBS unitary transformation and 
the four-photon coincidence postselection gives the two-ququart entangled state at the output in the form, 

|^ s = Q)4|o) c |o) d + Cl ci|i) c |i) d 

+ c 2 4|2) c |2) d + C 3C 3 |3) c |3) d 
+ 2(c o c 3 + Cl c 2 )\0) c \3) d 

+ 2(44 + 44)|3) c |o) d . 

The output state in eq. (|2"2")) , which is not that easy to manipulate experimentally, is the result for the most general 
superposition input state (pure state) with no coefficients equal to zero. With less general input ququart states, we 
can obtain a two-ququart entangled state which is experimentally more maneageable. 

Consider the case where c\ = c 2 = and c[ = 4 = in eq. ©■ Then the real SPDC input state to the PBS is 
given as eq. (J7J) with \ip) a = c \0) a + c 3 |3) a and \ip)b — c' \0)b + c' 3 \3)b- The postselected output state of the PBS then 
has the simple form 

\tp)out = co4|0)c|0)d + c 3 4|3) c |3) rf 

+ 2c c 3 |0) c |3) d + 244|3) c |0) d . 1 ' 



Note that, with the superposition of the ququart basis states |0) and |3) as the input, the output two-quauart entangled 
state also contains only these basis states. 
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C. Dichroic beam splitter 

We now consider two biphotons ququarts interfering at a dichroic beam splitter (DBS), see Fig. [2j The DBS is 
designed to reflect the Ai photon and to transmit A 2 photon, as in the case of the single biphoton ququart detection 
scheme in Fig. [TJ 

Depending on the frequency (Ai or A2) and the input spatial mode (a or b), the single photon entering the DBS is 
routed as follows, 

fljAi ^ dj\ l , a,jx 2 — > Cj\ 2 , 
bj\ 1 — > Cj\ 1 , bjx 2 — > rfjx 2 . 

Here the subscript j refer to the polarization state of a photon (if or V). 

Since the DBS always reflects Ai photon and transmits A2 photon regardless of the input spatial mode, Ai and A2 
photons behave differently at the DBS. For the Ai photon, the unitary transformation matrix is given as, 

10 



0100 



while for the A 2 photon, the corresponding unitary matrix is a different one: 

/ 1 

"2"- SJ! Si. < 26 > 

\Q 1 

These two matrices are simple mathematical representations of the wavelength dependent routing property of the 
DBS used in this experiment. 

The most general unitary transformation for the DBS that applies to the input four-photon or two-ququart state 
in eq. (O is therefore, 

jjDBS = jjDBS g jjDBS g jjDBS g jjDBS (2?) 

Again, this is a 256 x 256 matrix. 

For the general two-ququart pure state defined in eq. (J7J) as the input to the DBS, the output state can be determined 
by using eq. (j27|) and it is given as 

\^)out S = co^|0) c |0) d + C»,Ci|0)c|l>d + c o c' 2 |2)c|0) d + co4|2) c |l) d 
+ c lC ' a \l) c \0) d + Cl ci|l) c |l) d + Cl c 2 |3) c |0) d + ci4|3) c |l) d 
+ c 2 c[ ) |0) c |2) d + c 2 c / 1 |0) c |3) d + c 2 4|2) c |2) d + c 2 c / 3|2) c |3) d 
+ C3cf,|l) c |2) d + C3c£|l>c|3)d + c 3 c 2 |3) c |2) d + c 3 4|3) c |3) d . 

Equation (|28p has a few interesting properties. First, the 16 terms in eq. (f2"5f are all of the non-vanishing amplitudes. 
Therefore postselection is not necessary if only the two-photon event is present at the input, i.e., only the first term 
in eq. ([7]). Clearly, this is not possible with SPDC. However, the non-postselection nature of the DBS two-ququart 
entanglement scheme could prove to be important if recently reported semiconductor quantum-dot photon pair sources 
with supressed the multi-pair events become more widely available [3lJ. Second, the second and the third terms in 
eq. (|7J| do not result in any detectable two-ququart amplitudes in modes c and d so that they do not contribute to 
the final postselected amplitudes. With postselection, therefore, the SPDC input state could be considered as the 
ideal input as in eq. ^ in the case of DBS. Finally, if the input ququart states are formed by equal superposition of 
all ququart bases states, Co = c\ = c 2 = C3 and c' = c[ = c' 2 = d 3 , the output state turns out to be a two-ququart 
product state, i.e., \ip) out = |^)c ® \^)d- 

The DBS can be used to generate a number of interesting two-ququart entangled states starting from rather simple 
looking two-ququart product states as the input. For example, consider two input ququarts, 

\^ a = J-(\0) a + e^\3) a ), 

f (29) 

|^)6 = -^(10)6 + 6^13)6), 



(a) (b) 
V 



I t 



■ d' 

d 

FIG. 3: Scheme for linear optically entangling three (a) and four (b) biphoton ququarts. To generate the entangled state of n 
biphoton ququarts, n — 1 DBS would be needed. 

where tp± and (f2 are the relative phase terms. 

The DBS output state is a two-ququart entangled state of the form, 

|^)o U * = i(|0) c |0) d + e ^|l>c|2> d + e ^|2) c |l) d 



e 



|3)c|3) t 



IV. ENTANGLING SCHEMES FOR MULTIPLE QUQUARTS: EXAMPLES 

It is not difficult to see that the linear optical scheme for entangling two biphoton ququarts we have discussed so 
far can, in principle, be scaled up to entangle multiple biphoton ququarts. To discuss this possibility, we consider two 
simple examples: the three-ququart entanglement and the four-ququart entanglement. 

First, consider the three-ququart entanglement scheme shown in Fig. [3Ja). The three individual ququarts in the 
input modes of the DBS are 

Wa = |0) a , 

!</>)& = 10)6 + 13)6, (31) 
|^)c = |0)c+|3) c . 

Two individual ququarts in modes a and b get entangled by the action of the DBS and then one of the subsystem of 
the two-ququart entangled state becomes entangled with another ququart in mode c with a different DBS. 
As the result, the output state in modes a', 6', and c! is given as 

|^U = |0) a /|0)6'|0)c'+|0) a /|2) 6 Hl)c' 

+ |2) a ,|l)6<|0)c' + |2) a ,|3)6<|l) c >, [ ' 

which is clearly a three-ququart entangled state. 

The four-ququart entangled state can be generated in a similar way but with one more DBS. Roughly speaking, to 
entangle n biphoton ququarts, n — 1 DBS would be needed. Consider the scheme shown in Fig.[3jb) and the following 
four individual biphoton ququarts at the input modes of the DBS set: 



\i>)a = |0)a, 
|V)6 = |0)6 + |3) 6 , 
|V)c = |0) c +|3) c , 
= |0)«J. 

The output state is found to be 

\i>)out = |0)a'|0)6'|0) c '|0)^ + |0) a ,|l)6<|0) c '|2) d 



(33) 



|2) a H0)6'|l)c|0)^ + |2) a ,|l)6'|l) c H2)^, (34) 



which is a rather simple looking four-ququart entangled state. 
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To generate a more complex four-ququart entangled state, one simply needs to change one of the input states to a 
superposition state. For example, for the following biphoton ququarts entering the DBS set shown in Fig. EJb), 



the output state is calculated to be 



>« = |0)« + |3) a , 
)b = \0) b , 

>c = |0)c + |3) o , 

> c -|0> d + |3) d , 



)out = |0) a .|0) 6 '|0) c '|0)^ + |0) a -|0) b ,|2)c|l>d 
+ |l)a'|2)y|0) o ,|0) cl / + |l) a /|2) v |2) c /|l) tf 



(35) 



+ |2) Q ,|0MlM0) d ' + |2) a n0V|3) c ,|l)d' (36) 
+ |3) o /|2) y |l) cf |0) cl / + |3) a /|2) y |3)^|l) <l /. 

What we have shown in this section is that, with a relatively simple set of biphoton ququarts, it is possible to 
generate interesting-looking postselected entangled states of multiple ququarts. Scaling up to multiple ququarts 
becomes more or less a practical concern. 



V. DISCUSSION 



We have investigated scalable linear optical schemes to generate entanglement of multiple ququarts, in which the 
individual ququarts are constructed with the biphoton polarization states of frequency-nondegenerate SPDC. We have 
shown that it is possible to generate various entangled states of multiple ququarts by interfering them at beam splitter, 
such as, BS, PBS, and DBS. 

It is now relavant and interesting to discuss the efficiency of beamsplitter entangling scheme. Especially, we would 
like to know which kind of beamsplitter is best suited for generating an entangled state of multiple ququarts. 

Consider, first, the output two-ququart state of a BS shown in eq. (fl4)) . As discussed in section Ull A[ when the 
input state is of the most general form, there are total of 448 additional biphoton amplitudes that are omitted due 
to the two-ququart detection postselection. Experimentally, this means that there are simply too many amplitudes 
that do not result in four-fold coincidence so the experiment will be extremely inefficient in terms of the observable 
four-fold coincidence rate versus the SPDC pump power. 

The situation is much better with the PBS output state shown in eq. (|22p as there are only 36 biphoton amplitudes 
that cannot be expressed, hence not shown, in the ququart basis states defined in eq. ([T|). With the DBS, the double- 
pair terms in eq. ([7]) do not contribute any detectable amplitudes to the postselected output state in eq. (|28p. These 
omitted terms add up to total of 32 amplitudes and they cannot be detected by the two-ququart detection scheme. 

From the experimental point of view, it is best if there are no wasted amplitudes, i.e., the amplitudes that are thrown 
out by the four- fold coincidence circuit. This goal clearly is unachievable with SPDC sources but we may choose to 
improve the entangled state generation efficiency as much as possible by choosing to work with either the PBS or the 
DBS entangling scheme as both schemes offer roughly the same postselected-to- wasted amplitude ratio (1/3 for the 
PBS versus 1/2 for the DBS). One additional advantage of the DBS scheme would be that, with postselection, the 
input state could be considered as the ideal input as in eq. since the double-pair terms do not affect the output 
state. 

Finally, we note that, although the ideal input state of the form in eq. ^ is not achievable with SPDC, it may be 
possible to use recently reported semiconductor quantum-dot photon pair sources with supressed multi-pair events 
to prepare such an ideal input state [3l| . In this case, the non-postselection nature of the DBS entanglement scheme 
could prove to be important. 
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